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Summary. — The Standard Quantum Limit (SQL) for the measurement of a free 
mass position is illustrated, along with two necessary conditions for breaching it. 
A measurement scheme that overcomes the SQL is engineered. It can be achieved 
in three-steps: i) a pre-squeezing stage; ii) a standard von Neumann measurement 
with momentum-position object-probe interaction and lii) a feedback. Advantages 
and limitations of this scheme are discussed. It is shown that all of the three steps 
are needed in order to overcome the SQL. In particular, the von Neumann inter- 
action is crucial in getting the right state reduction, whereas other experimentally 
achievable Hamiltonians, as, for example, the radiation-pressure interaction, lead to 
state reductions that on the average cannot overcome the SQL. 

PACS 03.65. Bz, 42.50.Dv, 42.50.-p - . 



1. — Introduction 

The problem of achieving a sequence of measurements of the position of a free mass 
with arbitrary precision has received much attention in the past, especially as a tool 
for monitoring the presence of an external classical field weakly interacting with the 
mass itself — typically, gravitational waves [1, 2]. Originally, a Standard Quantum Limit 
(SQL) was devised [1, 3] stating that if a mass evolves freely for a time interval tf 
between two measurements, the uncertainty of the second measurement cannot he lower 
than ^"sQL ~ TUf/m, m being the inertial mass of the freely moving body. 

The SQL can be shortly illustrated as follows. Let us consider that the moving mass 
after the first measurement at t = is described by the following minimum uncertainty 
wave-packet (MUW) centered at q^, and moving to the right with momentum hk^ 
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If such wave-packet undergoes a free evolution for a time tf, its initial position variance 
(5^ = {Aq^{0)) increases as follows [4] 



where (• • •) = Tr[- • • g] denotes the ensemble average and AO = O — (O) for any operator 
O. Minimizing {Aq^{tf)) with respect to 5'^ in Eq. (2), one obtains {Aq^{tf)) > A^q]^. 
Thus, there is a limit on the accuracy of a subsequent position measurement, which 
originates from the spreading of the free mass wave-function. 

Is it possible to beat the SQL by preparing the moving mass in a different wave- 
function? The answer was given by Yuen [5], who identified a class of states — the con- 
tractive states — that have position uncertainty that decreases versus time, and thus can 
go below the SQL. Yuen also presented some measurement models of the Arthurs-Kelly 
type [6] that realize in different ways the state reduction toward contractive states [7, 8]. 
Ozawa subsequently proposed a different measurement scheme [9] that leaves the mov- 
ing mass in a contractive state, and can overcome the SQL. Moreover, after identifying a 
necessary condition for breaching the SQL [10, 11], Ozawa found a class of Hamiltonians 
that satisfy such condition [12]. 

In this paper we present a way of engineering ab initio a measurement scheme that 
beats the SQL. The scheme consists of three-steps: a prc-squeezing, a von Neumann 
measurement, and a feedback. It turns out that our measurement is equivalent (i. e. 
it has the same outcome probability distribution and the same state reduction) to a 
model belonging to a general class previously studied by Ozawa [12]. The object-probe 
interaction of the von Neumann measurement is Hj = qP. In our notation, lower-case 
operators denote system observables (the moving mass) and capital operators denote 
observablcs of the probe which, in our case, is a single mode of the electromagnetic 
field. Thus q denotes the position operator, whereas the role of the linear momentum 
P for the field is played by the quadrature P = = {A^e^'^ + Ae~'^'^)/2 of the probing 
mode, with annihilation and creation operators A and . The qP interaction may be 
difficult to achieve; hence we analyze also the case of a von Neumann measurement based 
on the interaction Hj = qA^ A, which is just the radiation-pressure Hamiltonian of an 
interferomctric measurement of a moving mirror position. However, in this case we show 
that we can only reach the SQL, but we cannot overcome it. 

The outline of our paper is as follows. Section 2 is a brief review of the formal 
framework for repeated quantum measurements. After giving the concepts of probabil- 
ity operator- valued measures (POM's) and instruments, we recall the notions oi precision 
and posterior deviation to describe the noise from the measurement device and the distur- 
bance from the state reduction. With these two concepts in mind we can recall the precise 
statement of the SQL due to Ozawa [10], who has given a necessary condition to breach 
the SQL. In Section 3 we see that the SQL can be overcome by a Gordon-Louisell (GL) 
measurement [13], as also shown by Yuen [8, 7]. After proving that every GL state reduc- 
tion can be obtained by means of a suitable feedback mechanism, we show how in this 
way it is possible to engineer a measurement scheme that beats the SQL. As announced, 
the scheme consists of the sequence of a pre-squeezing, a von Neumann measurement, 
and a feedback. In Section 4 we briefly recall the Ozawa's measurement models that 
satisfy Ozawa's condition, and show how our scheme realizes some measurements in this 
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class. In Section 5 we analyze the case of a von Neumann measurement achieved with 
the interaction Hi = qA^A. Section 6 closes the paper with some concluding remarks. 

2. — Repeated and approximated measurements 

In this Section we review the main points of the theory of repeated quantum mea- 
surements, and wc resume the formulation of the SQL given by Ozawa [10, 11], based on 
a necessary condition to breach the SQL. We will not give a complete general treatment 
of the subject, but only introduce the main concepts and notation that will be used in 
the following sections: for more extensive and rigorous treatments see Refs. [10, 11, 14] 
and references therein. 

2"1. POM's and Instruments. - A complete description of a quantum measurement 

consists of both: i) the probability density p{x\g)dx of the result x of the measurement 
when the quantum system is in the state described by the density matrix g; and ii) the 
state reduction Qx for the system immediately after the measurement with outcome 
X. Both p{x\Q)dx and Qx can be expressed in terms of one map, the so-called instrument 
dl{x), which is a linear map on the space of trace class operators g dl{x)g given by 

(3) p{x\g)dx = TT[dI{x)g] , 

dl{x)g 



(4) 



TT[dI{x)g] 



Hence, the instrument dl{x) gives a complete description of the quantum measurement. 
Sometimes, however, one is interested only in the probability density p{x\g) of the mea- 
sure outcome x, ignoring the state reduction: in this case it is sufficient to know the 
probability operator-valued measure (POM) dll{x) of the measurement, which provides 
the probability distribution of the readout x for any state g as follows 

(5) p{x\g)dx = TT:[gdn{x)] . 

By comparing Eq. (5) with Eq. (3) one can see that for every instrument dl{x) the 
corresponding POM dli(x) is defined through the trace-duality relation Tr[p(iII(a;)] = 
'Tr:[dl{x)g\. The correspondence between dl{x) and dXi{x) is not one-to-one, because the 
same POM can be achieved by different instruments dl{x), namely with different state 
reductions. 

Insofar we have given an abstract description of the quantum measurement, with no 

mention to the physical realization of the measuring apparatus. In order to have an 
output state that depends on the state before the measurement, the measurement appa- 
ratus must involve a probe that interacts with the system, and later is measured to yield 
information on the system. This indirect measurement scheme is completely specified 
once the following ingredients are given: i) the unitary operator U that describes the 
system-probe interaction; ii) the state |<^) of the probe before the interaction (we restrict 
our attention to the case of pure-state preparation of the probe) ; iii) the observable X 
which is measured on the probe. At the end of the system-probe interaction it is pos- 
sible to consider a subsequent measurement of a (generally different) observable Y on 
the system (with outcome y). Then, it can be easily shown (see, for example, Ref. [14]) 
that the conditional probability density p(y\x) of getting the result y from the second 
measurement — being x the result of the first one — can be written in terms of the Horn's 



4 



G. M. D'ARIANO, M. F. SACCHI AND R. SENO 



rule p{y\x)dy = {y\Qx\y) upon defining a "reduced state" Qx as in Eq. (4), where the 
instrument and the POM are given by 

(6) (m{x) = dxn\x)h{x) , 

(7) dI{x)Q = dxCL{x)Qh\x) , 

and the operator (l{x), which acts on the Hilbert space of the system only, is defined by 
the following matrix element on the probe Hilbert space 

(8) Cl{x) = {x\U\^) , 

\x) being the eigenvector of the observable X corresponding to eigenvalue x. As regards 
the evolution operator U, one can neglect, for simplicity, the free evolution during the 
measurement interaction time, and consider an impulsive interaction Hamiltonian that 
is switched on only for a very short time interval r with a very large coupling constant 
K, such that Kt is finite. For simplicity of notation, in the following we will implicitly 
include Kt in the definition itself of the interaction Hamiltonian. 

2'2. Precision and posterior deviation. - In any scheme for a quantum measurement, 

in principle there arc always two kinds of noise: i) the quantum noise of the observable 
that is intrinsic of the quantum state g, which is given by the variance (Ag^); ii) the noise 
due to the measuring apparatus, which is generally non-ideal, for example, because of a 
non-unit quantum efficiency, or as a consequence of the noise due to a joint measiircmcnt . 
This leads to an output probability distribution p{x\g)dx that is broader than {Aq'^). A 
POM c?II(a;) is said to be compatible with an observable q (or g-compatible) if it satisfies 
the relation 

(9) [dn(x),g]=0, 
namely the POM has the same spectral decomposition of q 

(10) dn{x) = G{x, q)dx = dx j G{x, q)\q){q\ dq , 

G{x, q) playing the role of a conditional probability density for the output x, given that 
the position of the system was q. 

The extrinsic instrumental noise or precision e^[Q] of the apparatus with g-compatible 
POM dll{x) for a measurement of the observable q, estimates the broadening of the 
intrinsic noise due to the measurement, and is defined as follows 

(11) e^[Q\ - 1 j{x- qfl^[QdR{x)\q){q\] dq = J Tr[{x - qfm{x)] . 



The overall noise or total uncertainty Aa;^ [g] of the measurement is the variance of the 
experimental probability distribution, namely 



(12) 
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where E[5(a;)||^] = / g{x)p{x\Q)dx denotes the experimental expectation value of the 
function g{x) of the random outcome x. The total uncertainty in Eq. (12) can be simply 
written as the sum of e^[g] and (Ag^) if the POM is unbiased, namely E[a;||^] = {q) for 
every state g. In that case one has 

(13) {q) = Tr[qQ] = jdxxp{x\g) , 

and hence q can be spectrally decomposed in terms of the POM itself, namely 

(14) q = J xdtl{x) ^ E[x\\g] = (q) , 



for all states g with Aa;2[^] < oo. The precision of a measurement realized with a 
g-compatible POM is equal to zero if and only if G{x,q) = 6{x — q), namely if the 
measurement apparatus is noiseless. Hence 

(15) e'^[g] =0 ^ dil{x) = dx\x){x\ . 

In this case the total uncertainty of the measurement is just the variance of the system 
state. 

For a g-compatible and unbiased POM, Eq. (12) can be rewritten as 



(16) Ax^g]=e'[g] + {Aq') , 

namely the intrinsic and the instrumental noises behave additively. The precision [g] of 
the measurement characterizes only the POM dll(x), namely the probability distribution 
of the measurement. On the other hand, the noise after the state reduction can be 
quantified by the posterior deviation [g] of the instrument /, which is defined as follows 



(17) 



'^^[0]= J j{q-xf{q\g^\q)p{x\g)dxdq = J Tv[{q - xf dl (x) g] . 



In Refs. [10] and [11] the quantity a^[g] is named "resolution": here we suggest to adopt 
the nomenclature "posterior deviation", as it is clearly a property of the reduced state 
gx- For g satisfying Aa;2[^] < 00 and Tr [A(f J dl{x)g] < 00 one has [11] 

(18) cj^[g] = j {Af)xp{x\g)dx + j [{q)^ - xf p{x\g)dx , 

where (• • ■)x = Tr[- • • p.,;] denotes the conditional expectation. Therefore for unbiased 
reduction — i.e. the average value of q after the state reduction is still equal to the 
outcome x of the measurement — cr'^[g] is just the variance of the reduced state averaged 
over all the readouts x, namely 

(19) Ti[qgx]=x =^ a^[g]= j{Af)xp{x\g)dx . 



6 



G. M. D'ARIANO, M. F. SACCHI AND R. SENO 



2'3. Generalized Standard Quantum Limit. - Let us consider a g-compatible and un- 
biased POM. At t = we measure the position of the system and we get outcome x and 
state reduction g^, with probability distribution Then wc let the system evolve 

freely for a time interval tf and perform a second measurement on Qx{tf)- According 
to a mean-value strategy, one predicts the second measurement to have the outcome 
h{x) = Tr[Qx{tf)q\ with an uncertainty given by 



(20) A{tf,g,x) = 



-,1/2 



x' — h{x)] p{x'\Qx{tf))dx' 



which, for an unbiased POM is just given by 



(21) A'{tf,Q,x) = Ax^Qx{tf)]. 

The "predictive uncertainty" A'^(tf,g) of the repeated measurement [10] for prior state 
Q is defined as the average of A^{tf, g, x) over all the outcomes x a,t t = 0, namely 

(22) ^^{tf,0) = j dxA^{tf,Q,x)p{x\g) . 

Ozawa [11] introduced a precise definition of the SQL as the lower bound of the predictive 
uncertainty A'^{tf,g) with the hypothesis that the averaged precision of the evolved 
state — the precision of a second measurement on the reduced state at time tf — is greater 
than the posterior deviation. Hence, for unbiased measurements satisfying the inequality 

(23) J dxp{x\g)e^[gx{tf)]>a^[g], 
Ozawa proved the bound 

(24) A'itf, g) > |Tr [^^[^(0), ^(i/)]] \' = A%q^ , 

where for the free- mass evolution Aqqj^ = htf/m. Eq. (24) is the general version of the 
SQL due to Ozawa. Using Eq. (21) one has 



(25) A^itf, g) = j dxp{x\g) e'[^.(t/)] + j dxp{x\g) TT[Aq^ gx{t f)]) > A 



2 

SQL ■ 



Therefore, a necessary condition for breaching the SQL corresponds to negating the 
hypothesis (23), namely that the precision of the measurement is less than the posterior 
deviation. We will call this condition ONC (Ozawa necessary condition). We recall that, 
while the precision (11) is a measure of the added noise and is related to the POM, the 
posterior deviation (17) is related to the state reduction. As an example, it is interesting 
to evaluate both cr^ [g\ and [g\ for a generalized von Neumann measurement model with 
interaction Hamiltonian of the form Hj = qO, where O is a generic observable of the 
probe, and q is just the quantity of the system that we want to measure — the position 
in our case (in the "standard" von Neumann model [15] O = P, where P is the linear 
momentum of the probe). From Eqs. (6), (7) and (8), it is apparent that both the POM 
and the instrument are only functions of the operator q, so that the precision (11) and 
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the posterior deviation (17) coincide. This means that condition (23) is verified only 
with the equal sign, and the von Neumann model just achieves the SQL. 

3. — Evolution operator for beating the SQL 

In this Section wc will resume the possibility to overcome the SQL by reducing the 
system to a contractive state. Then we will design an evolution operator IJ that realizes 
such reduction. Finally, we will generalize [/ to a class of operators which fulfill ONC. 

3'1. The contractive states. - The original argument for the validity of the SQL due to 
Braginskii and Caves [1, 3], was the following. Consider a mass m, whose position q has 
been measured once within a certain precision, and let it freely evolve for a time tf. In the 
Heisenberg picture the position operator evolves classically, as q{tf) = g(0) + p(jS)tf /m. 
From the Heisenberg uncertainty principle, one obtains the following position variance 
constraint 

(Ag2(<^.))^(Ag2(o)) + (Ap2(0)) {^-^^ 



However, as pointed out by Yuen [5, 8], this is not the correct general expression for 
{Af{tf)), because it neglects the correlation term 2Re(Ag(0)Ap(0)) = (Ag(O)Ap(O) + 
Ap(0)Ag(0)), hence implicitly assuming that it is greater than or equal to zero. The 
complete expression for the position uncertainty is 

(27) {Afitf)) = (Ag2(0)) + 2Re(Ag(0)Ap(0))^ + {Af{0)) (^^) . 

As a matter of fact, there is a class of states, the contractive states (CS), which have a 
negative correlation term in Eq. (27). When the reduced state for the free mass after 
the first measurement is a CS, {Aq^{tf))a: decreases in time before reaching a minimum 
value. If the measurement is sufficiently precise, the SQL can be overcome. In fact, 
if {Aq^{tf))x < (Ag^(0))x for every outcome x, then this is true also on the average, 
namely 

(28) J dxp{x\g) {Afitf)), < J dxp{x\Q) {Af{Q))^ . 
Eqs. (16) and (21) leads to 

(29) ^^{tf,Q) = j dxp{x\Q)e^[Q^{tf)] + j dxp{x\Q) {Af{tf))^ , 
which, for noiseless measurements e'^[Qx{tf)] = 0, gives 

(30) A\tf, q)= j dxp{x\Q) {Afitf)), < J dxp{x\g) {Af{0)), . 
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A particular example of a CS are the squeezed states, or twisted coherent state (TCS)[16]. 
A TCS \iJ,uaui) for the mass has the following position representation 

, . , / mw V r mwl + 2i^ n2 , « / ^" 

with fluctuations 

where — = 1, ^ = Im(/i*j/) > and a — + ipo, qo and po being real. The 
correlation function Re(Ag(0)Ap(0)) = —£,fi, is negative for ^ > 0, so that the position 
uncertainty {Aq^{tf)) decreases in time till, at tM = 2^/a-'|/i + z/p it reaches the minimum 
value (Aq^(tM)) — {A^)^^{htM /m) < A|q^ for sufficiently large ^ [10]. This implies that, 
for noiseless detection [Eq. (15)] the total uncertainty of the measurement goes below 
the SQL. A detection scheme of this type is described by the measurement with state 
reduction Q ^ Qq = \iivQ(jj) {ij,vQijo\ and probability distribution p{Q\q) = {Q\q\Q), Q 
being the output of a single measurement. This is a GL measurement with reduction 
operator d given by 

(33) n{Q) = \nuQu}){Q\ . 

It is easily proved that this scheme satisfies the request of unbiasedness for both the POM 
and the state reduction, and, as expected, Eq. (23) is contradicted, because cr^[^] = 
JdQp{Q\g)TT[AegQm > J dQp{Q\g)e^[gQ{tf)] ^ 0. 

Yuen [7, 8] and Ozawa have given different interaction Hamiltonians that realize the 
operator Q{Q) in Eq. (33). In the next Section, we will engineer ab initio a measurement 
scheme that has the state reduction \ij,uQu}). 

3'2. Gordon-Louisell measurements. - The definition of a general GL measurement is 
given by a state reduction operator of the form 

(34) h{x) = \^,){e,\ , 

where \0x) is a complete set, i. e. / dx\9x){0x\ = 1 and 1-02;) is a normalized physical 
state. The operator Q{x) in Eq. (34) abstractly represents an indirect measurement that 
leaves the system in the state \tpx), independently on the input state, and has probability 
distribution for the outcome x given by p{x\g) = {6x\g\0x)- In the particular case that 
\9x) is also an orthogonal set {Ox\Ox) = S{x — x'), there is a unitary operator U on the 
system-probe Hilbert space Hs ®'Hp that gives ^{x) in Eq. (34) through Eq. (8). In 
fact, let us define the self-adjoint operator 9 such that 9\9x) = 9x\9x)- Then, let us 
consider an indirect measurement scheme for 9, with the probe prepared in the (pure) 
state \(p). Comparing Eq. (8) with Eq. (34), we see that U must satisfy the identity 



(35) {9y\ {9x\U\9,) ® \ip) = {9y\i>x){0x\0,) = {9y\ilJx)S{x - z) 
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Hence, a unitary operator U that has the matrix elements (35) can be chosen as achieving 
the following linear transformation on Hs "Hp 



(36) 



How can we design such an operator tjl Consider the operator R that effects a rotation 
by 7r/2 in the {q, Q) plane, namely 



(37) 



R = exp 



2h 



ipQ - qP) 



Apart from an inversion and a trivial overall phase factor, the operator R corresponds 
to a mode-permutation operator since one has 



(38) 



^Vi(9)V'2((3)=V'2(-9)V'i(Q)- 



Then we introduce the notion of feedback operator F{x) G Tis, namely a self-adjoint 
operator that parametrically depends on the measure outcome x. This can be any self- 
adjoint operator fimction of x: however, we need to restrict the class of operator functions 
F{x) such that the integral Fs(X) = / dx F{x)®\9x){0x\ converges toawell defined self- 
adjoint operator acting on Hs ® "Hp, with X\9x) = x\9x) [17]. We choose the feedback 
Hamiltonian such that it connects the vectors |y) and |V'x)j namely 



(39) 



IV'x) = exp 



\v)- 



Finally, we write the operator U as follows 



(40) 



U = exp 



n 



FsiX) 



exp 



2h 



Now wc specialize Eq. (40) to the case described by Eq. (33). The wave- functions 
\ipx) and \0x) correspond to \jjLuQu;) and \Q) respectively, whereas the probe observable 9 
is the position Q, with outcome Q. This implies that the feedback operator has to shift 
the position of the system by a quantity equal to the output Q, namely Fs{Q) = pQ- 
Hence, the evolution operator has the form 



(41) 



U = exp 







n 


exp 



2h 



{pQ - qP) 



We emphasize that, as U contains a permutation operator, the initial probe state must be 
chosen to be a TCS. Hence, if the initial state of the system is we get R\il>)®\iivaijj) oc 
\^,v{—a)uj)®\'il)): in practice, the position of the object is indirectly squeezed by squeezing 
the probe position, and then exchanging the state of the system with that of the probe. 
For simplicity, we choose a = 0, which means that the initial position of the probe is equal 
to zero. After that, the feedback operator shifts the position of the system so that it finally 
corresponds to the output of the measurement, that is exp[—^F{Q)\\iivQuj) = \^vQuj). 
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Now we relax the value of the coupling constants in the interaction Hamiltonians and 
rewrite the evolution operator (41) in the general form 



(42) 



U = exp 



n 



exp 



ITT 

2h 



52 (s ^PQ - sqP) 



As regards the concrete feasibility of the evolution described by Eq. (42), the feedback 
part can be achieved by a position transducer, which displaces the mass by an amount 
giQ for every position outcome Q. Notice that wc do not need to realize the feedback 
Hamiltonian gipQ, because as a consequence of Eq. (8) this is equivalent to the system 
Hamiltonian gipQ, where the coupling is rescaled by the eigenvalue Q in place of the 
operator Q. The permutation operator will be discussed in the following section. 

4. — Comparison with previous models and realization 

In the first part of this Section we compare the unitary evolutions of the Ozawa's 
measurement models [12] satisfying ONC with our operator in Eq. (42). In the second 

part we factorizc our evolution into three steps, expressing the permutation operator 
by means of a pre-squeezing of both system and probe, followed by a von Neumann 
interaction Hj = qP. 

4"1. Ozawa's Hamiltonians and the conditions for beating the SQL. - In Ref. [12] 
Ozawa gives a class of measurement models that satisfy ONC. He explicitly writes the 
unitary system-probe evolution in an impulsive regime in terms of an interaction Hamil- 
tonian Hi, generalizing the standard von Neumann measurement model with Hj = qP 
[15]. The Hamiltonian is given by 



(43) 



Hi = -k+Qp - k-qP + k.iqp - QP) , 



with fc+ , and k^ as real parameters. For probe initial state with (fi{Q) = ip{—Q) 
and {(p\Q\(p) = (which is equivalent to condition a = in the previous section), and 
for the system initially in the state \tp), Eqs. (11) and (17) become 



(44) 
(45) 



e^[i;] = {l-cr{i;\f\i^)+d^{v\Q^\^) , 
a^[^] = {a - cfi^Plfli,) + {b- dfi^lQ'l^) . 



Here the coefRcicnts of the dynamics (a, b, c, d) — grouped in a matrix form that will be 
used in the following — are 



(46) 



a —b 




— c d 





coshlC + k^^^ 



sinh IC 



coshlC- k^^i^ 



where IC = y^k-^ + k+ fc_ can be either real or pure imaginary. 

The degrees of freedom of the above equations can be reduced by applying the criteria 
for a plausible object-probe interaction: one requests that if the uncertainty of the prior 
probe coordinate {(p\AQ'^\(fi) tends to zero, then the precision and the posterior deviation 
tend to zero. From Eqs. (44) and (45) it follows that one needs a = c = 1. It can be 
shown [12] that the condition c = 1 is also a consequence of the unbiasedness of the 
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POM. Moreover, as ad — 6c = 1, it turns out that b = d—1, and using Eqs. (46) one has 
d = 2 cosh /C — 1. With the above restrictions, the precision and the posterior deviation 

become = d'^ {(p\AQ'^\(p) and cr^ = {(p\AQ'^\(p) , independently on the initial system 
state. Thus, the ONC for beating the SQL can be expressed as \d\ < 1. In summary, we 
have 

a=l, b = d — 1 , c=l, 

(47) 

= |2cosh/C- 1| < 1 , ONC <^=^ < 1 . 

If /C is real, then d > 1; but, if /C is pure imaginary, d = 2 cos |/C| — 1, and —3 < d < 1. 
Thus, Ozawa's Hamiltonians satisfy ONC for certain values of the parameters {k+, fc_ , kz) 
which can be found by comparing Eqs. (47) with Eqs. (46), namely 

(48) k+=2kz, cosh/C = l + ^, -1< < . 

In particular, among the Hamiltonians (43), the model described by Ozawa in reference 

[9], for kz = 7r/3\/3 and k± = ±2^^, realizes the GL scheme \iii'Qll!){Q\ with d = 0. 
On the other hand, the original von Neumann model [15] corresponds to d = 1 and 

= CT^ = {(fi\AQ'^\ip). This means that there is a continuum of models for |d| < 1 that 
have a better precision than cr^ = {(p\AQ'^\ip) , and can possibly circumvent the SQL. 

We now look for the relation between our model and the model just described. We 
use the following realization of the angular momentum (complex) Lie algebra gl{2, C) 

(49) i^±(QP-m. 

One can easily verify that the above operators satisfy the gl{2, C) commutation relations 
[J+, J_] = 2Jz and [Jz, J±\ = ± j±, so that the exponential of their linear combinations 
can be faithfully represented by 2 x 2 Pauli matrices, independently of the value of the 
angular momentum J [the group gl{2, C) is complex, and the fact that the realization 
(49) does not preserve Hermitian conjugation is irrelevant for the group multiplication 
law]. Notice that the matrix in Eq. (46) is nothing but the Pauli representation of 
the evolution operator U = exp{—iHj /h) in terms of the gl{2,C) algebra realization 
(49). Both our operator U in Eq. (41) and Ozawa's in Eq. (46) for kz = 7r/3\/3 and 
k± = ±2kz, lead to the same coefficients (a, b, c, d), namely (1, —1, 1, 0). In fact, we have 
just seen that they both realize ^1{Q) in Eq. (33). In the general case, the relation 
between the coefficients of the operator (42) and the ones in Eq. (46) are 

a 

(50) 

c 

The last of Eqs. (50) clearly states that d belongs to the interval [—1,1], namely all 
possible models realizing the condition e < a are included. 



cos(— ) +5issm( — ) , b - 



1 . ,Trg2. , .7r52. 
■-sin(— )+5iC0s(— ) , 



ssin(— ), rf = cos(— ). 
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4'2. Realization through GL{2, C) elements. - In this Subsection we suggest a mea- 
surement scheme that reahzes the evolution operator in Eq. (42) or Eq. (46) through 
three steps in the GL{2, C) group. Wc arc describing an indirect measurement of the sys- 
tem position q through detection of the probe position Q, and the central element of our 
scheme is the von Neumann Hamiltonian J_ = iqP/h, which entangles the system object 
with the probe by shifting Q by q. Moreover, Eqs. (42) and (49), suggest that after any 
detection with output Q, a feedback of the form J_|_ = iQp/h is requested. This implies 
that, among the six evolution operators that can be obtained by permuting the expo- 
nentials of (J+, J_, Jz), we will not consider the three permutations with exp ^(^_|_J+^ 

applied before exp J_^: in the next Subsection we will analyze the remaining three 
cases. 

4'3. Feedback- assisted measurement. - Consider the sequence 

(51) exp (C+>^+) exp (C-'^-) exp (^2(^zJzj = exp (^A;+J+ -|- fc_J_ -|- 2kzJ^ . 

In this model, the operator exp {2C,zJz^ = exp iC,z{QP — QP)/^ does not entangle the 
system with the probe, but just pre-squeezes the states of both. Then, exp = 
exp {iC,-qP/?)}j entangles the probe with the system. Finally, the operator exp = 

exp (^i(^Qp/hj corresponds to a feedback mechanism that shifts the system position by 

C+Q. By means of the Lie algebra decomposition formulas [18], the relations between 
the parameters <^+, C-, a^^d the Ozawa's write as follows 



(52) 



(l + C+C-)ef% -& = C+e- 



C-e' 



d: 



The range of the coefficients {(+,(-, in which the SQL can be overcome is [see Eq. 
(47)] 



(53) 



C- 



G > 



Now a problem arises, regarding the class of Ozawa's Hamiltonians represented by this 
model. In fact, as d= e~^' > 0, it is clear that a part of the evolutions that circumvent 
the SQL are excluded by this parameterization. In particular, rewriting {k+,k-,kz) in 
terms of (C+, C-,Cz) as in (53), one gets 



=-G _ 1 



e-C- + 3 



1/2 



arcsm 



1/2- 



(54) 



k- = 



-Q-l 



1 k, ^ 

2 k- 



where a narrower interval for kz/k- than in Eq. (48) is obtained. One is lead to think 
that there must be some other schemes that give the Ozawa's Hamiltonian for d > Q. 

However, an outlook at the remaining two permutations with exp iC,_j\ acting before 
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exp reveals that also in these cases d — e > 0, which implies that a substantial 

part of the Ozawa's Hamiltonians cannot be realized with our scheme. Notice, however, 
that the Hamiltonian with coefficients = 7r/3\/3 , A:+ = 2kz , fc_ = —^kz, which beats 
the SQL, is included in our scheme (51), and can be achieved in the limit of infinite 
that is for very high squeezing for both probe and object preparations. The realization 
of the squeezing of a mass position remains as a challenge for experimentalists. We stress 
the fact that, if either (z or in Eq. (51) are set equal to zero and the conditions for 
breaching the SQL are imposed, the remaining two coefficients become equal to zero, 
too. That is, both the feedback and the dilatation are essential in order to achieve a 
detection scheme that possibly beats the SQL. 

Now we calculate the operator fl{Q) relative to the evolution operator U just de- 
scribed. Because of the peculiar form of the operator U, also the operator (l{Q) can be 
factorized into three parts. By separating the exponentials that act only on the system 
from those acting also on the probe, we obtain 

h{Q) = exp {iC+Qp/n){Q\ exp (^iC^qP/hj exp {iQ^QP/n) \ip) 

(55) X exp {-Kzqp/K) ■ 

By imposing the conditions (53) for a plausible measurement, the operator (55) becomes 

(56) n{Q) = exp [i (e«- - l) Qp/U] exp {-i(zqp/h)if[e'^^ {Q - q)] . 

In particular, for probe initial state \(p) chosen as a TCS \(p) = IfiuOu}), we obtain 

^{Q) = exp [i (e^" - l) Qp/h] exp (-iCzqp/h) 
57 \— ^ exp - — - ^e^- Q-g 2 



Thus, the POM becomes 



(58) 



dXKff) 



mw 



7rh\iJ, — p\' 



exp 



mw 1 

h \fi — u\ 



dQ , 



which is Gaussian and unbiased. Hence, the precision [g\ , is independent on the state q, 

and is given by the variance of the Gaussian, namely e^lg] = e'^G (^(p\AQ'^\ip). This is pre- 
cisely what we expected, because, as we have seen in Subsection 41, e'^lg] = d'^{ip\AQ'^\ip}. 
We recall that, for very high squeezing, the operator (57) reduces to the GL operator 
(33), and e^[g\ tends to zero. 



5. — Radiation-mirror interaction 



In our knowledge there is no viable way to achieve the von Neumann interaction 
H oc qP. In order to approach its behavior, other Hamiltonians have been suggested. 
In particular. Walls et al. [19] described the interaction Hi = Aq, where a harmonic 
oscillator, the end mirror of a cavity, is interacting with the radiation in the cavity 
through radiation pressure. In this Section, we examine the situation in which the mirror 
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can be considered as a free mass, similarly to the case of a typical gravitational waves 
interferometer, where the wave detector that is attached to the mirror is a very massive 
bar (m ^ lOOfcg [3]), which has negligible oscillation frequency. We rewrite the radiation 
interaction Hamiltonian in the form 

(59) Hi = -nKmqA^A . 

For the mirror at one hand of a cavity of length L, the coupling constant Km can be 
derived [19] as Km — ujq/L, where uio is the resonance frequency of the cavity. With the 
impulsive approximation K^t ~ 1 the unitary evolution operator becomes simply 

(60) U = exp{-iHiT/h) = expiiqA^A) , 

where the mass position q = q/l-r is rcscaled by Ir = L/uJ{)T. We know that the von 
Neumann-type measurement schemes allow only to reach — not to beat — the SQL, and 
from the previous sections we learned that in order to overcome the SQL we need to 
add a feedback and a prc-squcczing to the present scheme. Therefore, we propose the 
following simple detection scheme. The probe is prepared in a highly excited coherent 
state \ip) = \a), where a = —i\a\, and then we detect (by homodyning it) the scaled 
quadrature -^/|q!| = {A^ + A)/2\a\. The operator Cl in Eq. (8) becomes 



(61) fig(X) = ( ^ j exp 



-|a|2(X + ie'i)2-H!(i + e2'i) 



and for the POM we have 

(62) dtlgiX) = dX (^^^) ^ exp [-2\af{X - sinqf] . 

For states with (g) <c 1 and (Aq^) <C 1, such that sin^ can be approximated as sing ~ q, 
one has the unbiased POM 

(63) dtl^iX) = dX {^^^ ' exp [-2|a|2(X - ^^j . 

Hence, for small mirror displacements, we can check if there are violations of the SQL. 
Notice that, within the small q approximation, the operator (61) becomes 

(91 12 \ 1/4 
exp[-|a|2(X-g)2+i|a|2(g2x + g-2X)] . 

The feedback shifts the object position by a quantity proportional to the output X. The 
operators Oq(X) becomes 

(65) fi(X) = exp[i(T7r)X£] Vl^iX) , 

where p = p/pr, with pr = fil^^. As the feedback is represented by a unitary operator, 
it does not modify the POM, but changes only the state reduction, so that the posterior 
deviation (t[qY changes and condition (23) can be beaten. 
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Now we check if there is violation of the SQL. From Eq. (63), the precision of the 
measurement can be calculated as e\g]'^ = Z^/(4|ap), which is independent on the state 
of the system: however, the position fluctuation depends on it. Thus, we must choose 
the object initial wave-function at time t = —t before the interaction, then calculate 
the reduced state after the first measurement at t = 0, and then evolve it for a time 
interval tf. We choose tjj{q, — r) to be a MUW as in Eq. (1). The calculation can be 
made in two main steps. We first multiply the function flq{X) by ^j(q, — t) (the reduction 
operator Clq{X) depends only on the position operator, and iplq, — r) is written in the 
position representation). The second step consists in applying the feedback operator 
exp[i{T' /t)Xp\ to the result, shifting the g-coordinate by Irir' /t)X. The reduced wave- 
function 'ip{q, 0) turns out to be of the form 



(66) tp{q, 0) oc exp —q 

|2 



+2q 



\a\ 

Ir 



X 



a 



go ikp 
45g + 2 



L r 



X{l-iX)- 



\a\ t' IrX 



By remembering that the probe is prepared in a highly excited coherent state, and by 
keeping only the higher order terms in \a\, we get 



(67) 



4>(q, 0) oc exp ■ 



\a\^ 



{1 - iX) { q + -\a\lrX 



The comparison between Eq. (67) and Eqs. (31) and (32) shows that 'tp{q, 0) can be a TCS 
or not, depending on the value of the output X of the measurement. In fact, the variance 
(V'|Ag2(0)|i/)) = l'^/{4\a\'^) can be set equal to hlfx-ul"^ /{2mco), with \fj,\'^ - = 1, if \a\ 
and Ir are appropriately chosen. The time evolution for {ip\Aq'^ {t)\tp) is described by Eq. 
(32) with ^muj/iniiJ, - v]'^) = -\a\'^X/l1, which means that £, = Im(/i*i/) = -X/2. This 
implies that the reduced state (67) is a TCS only when the measurement result gives 
a negative value X < 0. However, this occurs only with 50% probability, because the 
probability density for X is Gaussian and centered in X = 0. Hence, on the average — i. 
e. in Eq. (29) — the SQL is not beaten. 

Even a pre-squeezing of the free mass before the measurement does not solve the 
matter in the hand. In fact, the operator Q,{X) in (65) would become 

(68) diX) = exp[i{T' /t)Xp] ClqiX) exp[-i{T" /T)iqp + pq)] . 

The squeezing operator eKp[~i{T" /T){qp+pq)] squeezes the g-coordinate as g ^ qe'^'^ 
This reflects on the POM (illq(X), which changes to (ill^gx"/x {X) so that the precision be- 
comes e^[^] = i^/(4|ape^'^ Z"^). This precision becomes very small both if the initial state 
for the probe is very excited and if the pre-squeezing on the system is very high. However, 
the squeezing operator does not modify the functional form of the initial MUW, because 
it only rescales its average values: in fact exp[— i(r"/r)(gp -I- pg)]|(7) = ^'^'^\e^ ^^g), 
which implies that ilj{q,—T) —^ il){qe^" /'^ ,—t). Thus, the reduced state of the system 
after the pre-squeezing, the measurement and the feedback is again of the form (67), 
and the SQL is not overcome. This result suggests some considerations about the way 
of beating the SQL. We know that, for generalized von Neumann measurements, the 
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averaged precision equals the posterior deviation, which, if the condition (19) is satisfied, 
coincides with the variance of the reduced wave-function averaged on the previous read- 
out. Thus, when starting from von Neumann schemes, in order to beat the SQL, wc must 
change (Ag^(O)). We accomphshed that by applying a pre-squeezing and a feedback to 
the system, and we have seen that, if the probe is initially in a coherent state, at best 
the SQL can be reached. The reason is that, oven though the precision tends to zero for 
very high squeezing, the uncertainty of the reduced state does not decrease with time, in 
general. This suggests us to squeeze also the probe ([9, 19]) as we have done in Section 
(3). This should reflect on the state of the free mass by reducing it to a contractive state, 
independently on the output of the measurement. Thus, we now calculate the operator fl 
when the probe is prepared in a squeezed coherent state |— z|a|,— r) = D{—i\a\)S{—r)\0). 
The reason why we do not take a squeezed vacuum as in Ref. [19] is apparent from the 
form of the evolution operator (60): it represents a rotation in the complex phase space, 
and if applied to the vacuum it does not modify the state. The evaluation of the operator 
dq{X) needs a quite long derivation: we only report the result 



(69) ng{X) = 



2^ 



exp < — ijap cosg(2X — sing) — |ap(X — sinqY 



tang 



1 



o2r 



1 



where 



(70) 
(71) 



-2r , 



K = e sm g -|- e cos q , 
4> = —arctan (e^'^tang) 



If we suppose that the mirror displacements arc very small, which means that (q) <C It, 
we have the first-order approximation for flq{X) 



(72) 



%(X) 
X exp 



2ldV'-\3 



{Xf+q-2X) - \afe^^ {X - q)' + %e^^ - 1) 



As we have seen in the previous derivation, the pre-squeezing does not change the form 
of the ^-representation of the system wave-function, but only rescales both the average 
and the variance by squeezing factors. Hence, wc choose an interaction described only 
by a measurement followed by a feedback, like in Eq. (65). Once again, we first multiply 
the initial wave-function by dg{X), then shift q by \a\lrX/T and finally impose that la] 
is very high. We get 

(73) V'(5,0)(xexp|-^P(^l-iJ) (^q + l^\a\lrX^ | , 

which is of the same functional form of the state reduction (67), and is contractive only 
for negative values of the output X of the measurement, that still occurs with 50% 
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probability. Therefore, by replacing the standard von Neumann Hamiltonian Hi = qP 
with the interaction of the free mass with the radiation pressure, we can beat the SQL 
for certain values of the output, but on the average wc only get a very narrow final 
state. This holds true even for a probe prepared in a highly squeezed state, and using a 
pre-squeezing and a feedback. 

6. — Conclusions 

In this paper, wc have engineered ah initio a measurement scheme that allows to beat 
the SQL. We have shown that our scheme belongs to a class of measurement models pre- 
viously studied by Ozawa. The measurement can be performed in three-steps, involving 
a prc-squcczing stage, a von Neumann interaction Hi = qP, and a feedback. For large 
squeezing, and with the probe prepared in a TCS, the state reduction puts the system 
into a TCS too, and the SQL is breached. When the standard von Neumann interaction 
is replaced by a mirror-radiation interaction Hj = qA'^A, the SQL cannot be beaten, 
even though all other conditions are kept the same. This is due to the fact that in the 
limit of small free mass displacements the radiation-pressure Hamiltonian gives the same 
probability distribution than the von Nciimann one, but not the same state reduction. 
In fact, the von Neumann interaction has the capability of transferring the shape of the 
wave-function from the probe to the system, so that the reduced system state can be 
made contractive, and thus the SQL is breached. On the other hand, this is no longer 
true for the radiation-pressure Hamiltonian, where the state reduction of a MUW is a 
very narrow state, which can be contractive with 50% probability, but on the average 
the SQL is not beaten. Thus we conclude that the experimental realization of the pre- 
cise form of the von Neumann Hamiltonian is in order, if the measurement apparatus is 
designed to beat the SQL. 
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